Near-degeneracy of several pairing channels in multiorbital models for the Fe-pnictides 
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Weak-coupling approaches to the pairing problem in the iron pnictide superconductors have predicted a wide 
variety of superconducting ground states. We argue here that this is due both to the inadequacy of certain ap- 
proximations to the effective low-energy band structure, and to the natural near-degeneracy of different pairing 
channels in superconductors with many distinct Fermi surface sheets. In particular, we review attempts to con- 
struct two-orbital effective band models, the argument for their fundamental inconsistency with the symmetry 
of these materials, and compare the dynamical susceptibilities of two and five-orbital tight-binding models. We 
then present results for the magnetic properties, pairing interactions, and pairing instabilities within a five-orbital 
tight-binding Random Phase Approximation model. We discuss the robustness of these results for different dop- 
ings, interaction strengths, and variations in band structure. Within the parameter space explored, an anisotropic, 
sign-changing s-wave (Aig) state and a (1^2 _y2 (Big) state are nearly degenerate, due to the near nesting of 
Fermi surface sheets. 

PACS numbers: 74.70.-b,74.25.Ha,74.25.Jb,74.25.Kc 



I. INTRODUCTION 

The undoped Fe-pnictides are semi-metallic materials 
which exhibit structural and spin density wave (SDW) anti- 
ferromagnetic transitions. Photoemission^ and density func- 
tional theory (DFT) calculations^'^ "^'^ find that over an energy 
range near the Fermi energy the electron bands are made of 
states from the Fe-pnictide layers (Fig. [T]) of predominantly 
Fe character. The two dimensional nature of these layers is 
such that the Fermi surfaces consist of hole cylinders around 
the r point and electron cylinders around the M point of 
the 2 Fe/cell Brillouin zone (see Fig. [2]). In the undoped 
system, the near nesting of the hole and electron Fermi sur- 
faces can give rise to a colinear antiferromagnetically ordered 
state within a weak coupling approximation^, and this state 
has been confirmed in neutron experiments^. In both the 
electron and hole doped cases, superconductivity appears in 
proximity to or coexists with the antiferromagnetic SDW or- 
^gjj7j8j9jio therefore natural to consider the possibility that 
spin fluctuations provide the pairing mechanism and various 
random phase^^'^^'^^'^"^'^^ (RPA), fluctuating exchangai^i^iS 
(FLEX) and renormalization group^^ calculations have been 
reported. These have made use of different approximations 
to the band structure and obtained a variety of different gap 
structures. 

At the same time, various experimental probes of differ- 
ent Fe-pnictides appear to indicate quite different symmetries 
of the order paramete r^^i^^i^^i^'^i^^i^^i^^i^^i^^i^^i^^i^^ ; as discussed 
below, there is a real possibility that different symmetries are 
realized in different materials. It is therefore important to in- 
vestigate the origin of the discrepancies among the various 
theories with regard to their predictions for the superconduct- 
ing ground state. In the case of the cuprates, where a single 
Cu-0 hybridized band predominates at the Fermi surface, dif- 
ferent methods and band structures all lead to (i-wave pairing. 
If the situation is qualitatively different in the Fe-pnictides, we 
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FIG. 1 : (Color online) The crystal structure of LaOFeAs showing the 
Fe-As layers with an Fe square lattice (red) and As atoms (yellow) in 
a pyramidal configuration above and below the Fe plane. 



need to isolate which of the apparently small differences in the 
various approaches is responsible for the different results ob- 
tained. The reward for such an effort may be considerable, 
since the evidence suggests that correlation effects in these 
materials may be modest compared to the cuprates, raising 
the possibility of a quantitative theory of the superconducting 
state. 

Here we describe the results of an RPA calculation of the 
pair coupling strength and the momentum dependence of the 
gap function for a multiorbital tight binding parametrization 
of a DFT band structure for the LaOFeAs material by Cao et 



FIG. 2: (Color online) The Fermi surface for LaOFFeAs found in 
Ref.i. 
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alA. We begin in SecHIlby discussing the two orbital model 
for the electronic structure of the pnictides that was introduced 
by some of the present authors and has been studied by var- 
ious groups, along with some of the inconsistencies associ- 
ated with this approach. In Sec. Hill we discuss our results for 
the electronic structure using a five-orbital fit to DPT. Then in 
Sec. [iVlwe examine the spin and orbital susceptibilities and 
the effective pairing interaction for both models. In Sec.|Vl an 
effective pairing interaction strength A and gap function g{k) 
are introduced and in Sec.|Vllwe discuss our results for A and 
g{k). In Sec. lVIIl we examine the spatial and orbital structure 
of the pairing interaction and the resulting pairs. Our conclu- 
sions and comparison with previous work are contained in the 
final section rvini 



II. TWO-ORBITAL MODEL 

Effective models of the band structure based entirely on 
Fe orbitals should be possible because DFT tells us that the 
states due to the As 4p orbitals are located approximately 2 
eV below the Fermi level. The As orbitals allow for hybridiza- 
tion with the Fe 3d states, however, and therefore an effective 
Fe-Fe hopping Hamiltonian can be constructed, provided the 
symmetries of the entire FeAs layer are respected. It is con- 
ceptually simplest to work therefore with a square lattice with 
sites corresponding to Fe atoms and introduce a set of effec- 
tive hoppings between these sites. If the primitive unit cell is 
taken to be a square containing a single Fe atom, the effec- 
tive Brillouin zone is the square shown in Fig. [3^. Since the 
true primitive unit cell contains two Fe atoms, the true Bril- 
louin zone is smaller by a factor of two, as shown in Fig. [Sj). 
The Fermi surface in the correct small zone should be ob- 
tained from the effective large zone by folding the effective 
zone about the black dashed lines in Fig. [3^. Sheets around 
the M point in Fig. [3^ thus correspond to sheets around the V 
point in Fig.[3j). 

In Qi et al}^, a 2-orbital tight binding model was used to 
carry out an RPA calculation of the spin and orbital fluctua- 
tion pairing interaction. Here we briefly describe this model 
and discuss why it is insufficient to approximate the full band 
structure, especially with regard to the correct orbital weights 
along the Fermi surface sheets. The Hamiltonian for the two 



FIG. 3: (Color online) (a) Effective Brillouin zone for Fe-Fe tight- 
binding model; (b) Correct zone. In (a) the dashed black line shows 
the boundaries of the correct zone, the dashed blue lines show the 
two disconnected pieces of the path that have to be added to give a 
corresponding path in (b). 

band model can be written as 

i^O = ^ ^ta [(e+(^) - 1 + e-(^)^3 + e:j,y{k)Ti\ ^Jka 
ka 

(1) 

Here, are Pauli matrices and i/j^^ = [d^^{k), d'^^{k)) is a 
two-component field which describes the two ''dxz" and ''dyz' 
orbitals. The energies e±{k) and exy{k) are parameterized in 
terms of four hopping parameters ti : 

e+(/c) = — {ti -\- 12) {cos kx -\- cos ky) — 4:ts cos kx cos ky 
e-{k) = —{t I — t2) {cos kx — cos ky) (2) 
^xy {k) = —4^4 sin kx sin ky 

Diagonalizing Hq one has 

^sa{k) = Y,at{k)j,^{k) (3) 



Here, jua{k) destroys an electron in tht u = ± band with 
spin a and 




With the addition of a chemical potential the Hamiltonian 
takes the following form 

^0 = E E (^-(^) - 7+ (fc)7..(fc) (5) 

ka v=± 

where the band energies are E±{k) = ^-\-{k) ± 
^ e^_{k) + ^ly{k). Measuring energy in units of for 
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ti = —1, t2 = 1.3, ts = = —0.85 and /i = 1.5 we 
exhibit the Fermi surface in Fig. (H displayed in the large ef- 
fective zone. To obtain the corresponding Fermi surface in the 
small zone, a folding across the line joining (tt, 0) and (0, tt) 
is required, as discussed above. 

The Fermi surfaces ai around (0, 0) and a2 around (tt, tt) 
are hole pockets associated with E- (kf) = and the f3i and 
(32 Fermi surfaces around (tt, 0) and (0, tt) are electron pock- 
ets from Ej^(kf) = 0. We note that the displacement of the 
a2 Fermi surface from the F-point to (tt, tt) is an artifact of 
the 2-orbital approximation. As shown by various DFT calcu- 
lations^'"^'^ and noted by Lee and Wen^^, the orbital states that 
have significant weight near the (3 Fermi surfaces include, in 
addition to the dxz and dyz orbitals, a dxy orbital. In addition, 
while the Fermi surface sheets shown in Fig. |4] fold down to 
give two hole Fermi surfaces around the F point of the 2Fe/cell 
Brillouin zone, there should in fact be two hole Fermi surfaces 
around the F point of the large, effective Brillouin zone. This 
is known from the wave functions found in the band structure 
calculations. Finally, the 2-orbital model lacks the flexibility 
to fit the Fermi velocities found in the band structure calcula- 
tions, giving Fermi velocities on the electron (3 Fermi surfaces 
that are anomalously small compared to the velocities on the 
hole Fermi surfaces. 



III. 5-ORBITAL MODEL 

In principle, one could capture the correct behavior near the 
a and (3 Fermi surface sheets by treating a 3-orbital d^z, dyz^ 
dxy model. However, with short range hoppings, this leads 
to the appearance of an extra unphysical Fermi surface and a 
fourth orbital is required to remove it^^. In addition, recent 
theoretical calculations using 5 Wannier d-orbits per Fe site 
find that one can obtain an excellent representation of the elec- 
tronic structure within a ±2eV window of the Fermi energy 
Furthermore, the values of the Coulomb interaction parame- 
ters obtained for the Wannier basis are such that the average 
Coulomb interaction is small compared to the bandwidth, im- 
plying that one is dealing with a weakly coupled system^^. 

At this point it therefore seems best to use all five d orbitals 
in developing a tight-binding modelii. Here, using a Slater- 
Koster based parametrization which respects the symmetry of 
the FeAs layers, we fit a five-orbital (dxz, dyz, dxy, dx2-y2, 
and dsz-^-r"^) tight binding model to the DFT band structure 
determined by Cao et al^. We will use an orbital basis that is 
aligned parallel to the nearest neighbor Fe-Fe direction rather 
than the Fe-As direction. With this choice we avoid the ne- 
cessity of introducing a second, rotated coordinate system in 
addition to the one that is used to describe the single Fe unit 
ceU. 

The Hamiltonian for the 5 band model takes the following 
form 

^m^mn) dl^^(k)dncj(k) (6) 

kcT mn 

Here d\^ ^{k) creates a particle with momentum k and spin 
(J in the orbital m. The kinetic energy terms (,mn{k) together 
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FIG. 4: (Color online) Fermi surface for the 2-orbital dxz-dyz model 
with /i = 1.5 showing the 0:2 FS sheet around k = (tv^tv) rather 
than around k = (0, 0). The orbital contributions of the band states 
that lie on the different FS sheets are shown color coded: dxz (red) 
and dyz (green). 



with the parameters for a 5 band tight binding fit of the DFT 
band structure by Cao et ah are listed in the appendix. A 
diagonalization of this Hamiltonian yields the eigenenergies 
and the matrix elements analogously to the two band case dis- 
cussed above. 

In Fig. [5] (a), we have plotted the resulting band structure 
in the backfolded "small" Brillouin zone while the Fermi sur- 
face sheets for zero doping are shown in Fig.[5](b). The col- 
ors correspond to the dominant orbital weight of each band 
in momentum space. The gray lines represent the DFT band 
structure by Cao et al and the comparison shows, that the 
5 band fit approximately reproduces the DFT bands, espe- 
cially in the vicinity of the Fermi level. It is obvious that the 
dxy contribution plays an important role in building up the 
electron-like Fermi surfaces {(3 sheets) around the M point 
of the "small" Brillouin zone. In the unfolded Brillouin zone 
(Fig.[5](b) the dxy orbital and dyz (dxz) orbital contribute the 
dominant weights to the band states at the (3i{(32) Fermi sur- 
face. To confirm this we also show the orbital weights as a 
function of the winding angle on the different Fermi surfaces 
in Fig.[6l 



IV. THE SPIN AND CHARGE SUSCEPTIBILITIES 

A. Noninteracting susceptibilities 

In this section we examine the RPA enhanced spin and 
charge susceptibilities for the multiorbital model that we in- 
troduced in the previous section. The spin operator for an 
orbital s is defined as 

Ss i<l) = \Yl + l)^-l^d,0 (k) (7) 

where a and P are spin indices. The spin susceptibility can 
then be calculated from the Matsubara spin- spin correlation 
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FIG. 6: (Color online) The orbital weights as a function of the wind- 
ing angle (j) on the different Fermi surface sheets. The different colors 
refer to dxz (red), dyz (green), dxy (yellow), dx2_y2 (blue), d32 2_^2 
(magenta). 
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FIG. 5: (Color online) (a) The backfolded band structure for the 5 
band model with F, X, and M denoting the symmetry points in the 
real Brillouin zone corresponding to the 2 Fe unit cell. The main 
orbital contributions are shown by the following colors: dxz (red), 
dyz (green), dxy (yellow), dx2_y2 (blue), d^z'^-r'^ (magenta), and 
a strongly hybridized dxz-dyz band (brown). The gray lines show 
the correct DFT band structure calculated by Cao et al (b) The FS 
sheets of the 5 band model for the undoped compound (x = 0). 



function 



{XiYdqM = 3 dT e'^- {TrSs{q,T)St{-q,Q)) (8) 

with r the imaginary time and uo a Matsubara frequency. In the 
same way we can define the Fourier component of the charge 
density for the orbital s as 



k,af3 

and we can calculate the charge susceptibility from 



(9) 



(X 



Jo 



dre'^^{Trns{q,r)nt{-q,0)) (10) 



In a more general formulation the susceptibilities are func- 
tions of four orbital indices. For the non-interacting case 



{Xo)st {Xi)st equivalent and can be written 



xinq 



= /"%re--^^(T,4Jfc,T) (11) 

kk' a(3 



dqa{k + q, r)dl^{k\ ^)dtp(k' - q, 0)) 



Now we can derive an explicit expression for the non- 
interacting susceptibilities from 



k,iuJn 



(12) 

where N is the number of Fe lattice sites, /3 = 1/T, and the 
spectral representation of the Green's function is given as 



Gsp{k,iuJn) = ^ 



<(fe)ar(fe) 



(13) 



Here the matrix elements a^^{k) = {s\fik) connect the orbital 
and the band space and are the components of the eigenvectors 
resulting from the diagonalization of the initial Hamiltonian. 
Performing the Matsubara frequency summation and setting 
iujn ^ uj -\- iO~^ we find the retarded susceptibility 



-y 



aUk)al*{k)al{k + q)al*{k + q) 



iV^ uj + E^{k + q)- E^{k)+iQ+ 



x[f{E,{k + q))-f{E^{k))] 



(14) 
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B. Random phase approximation for multiorbital system 



m 



(b) 



Now we consider Coulomb interactions of the electrons on 
the same Fe atom in an RPA framework. We distinguish be- 
tween an intraorbital interaction U of electrons in the same 
orbital and an interorbital interaction of electrons in different 
orbitals V. We can also take the Hund's rule coupling into 
account that favors the parallel alignment of electron spins on 
the same ion and is described by an energy J > and the pair 
hopping energy denoted by J^ These interactions are gener- 
ated automatically in multiorbital models with general two- 
body interactions using a Hubbard-type approach restricted to 
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. Thus in general one can write 



Hi, 



E 



E 



' 2 



E 



E 



where riis = rii^s] + ^isi- We have separated the intraor- 
bital exchange J and "pair hopping" term J' for generality, 
but note that if they are generated from a single two-body term 
they are related by J' = J/2. If we now split the interaction 
Hamiltonian into singlet (H-^]), triplet (H-^^) and pair (H^^l) 
channels, we find 



E 



V ,^ , ; 



and 



it) 



where a. 



= -E 

i,st 

2 Sis, as well as 



(16) 



(17) 



J' 



Htl = E E - ^MsAacu^cu. (18) 
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FIG. 7: (Color online) The bare spin susceptibility xs{q) calculated 
for the undoped 5 orbital model without (a,c) and with (b,d) the ma- 
trix elements. While (a) and (b) show cuts of the susceptibility along 
the main symmetry directions, (c) and (d) show the susceptibility in 
the first quadrant of the effective Brillouin zone. 



where a ^ b. Our notation of the interaction parameters U, 
V, J and J' can be compared to the notation in Kubo^^ as 



J' 



(21) 



where U, V, J and J' denote the interaction parameters in- 
troduced in Ref. [37[ In Fig. [7] we show results for the static, 
homogeneous bare spin susceptibility 



Xs 



(q,o) 



(22) 



The RR\ susceptibilities are obtained in the form of Dyson- 
type equations as 



IaO )st — Xst IaO )uv\^ )wzX. 



RPA\pqi 



wz 
Cst 



and 



RPA\pq, 



(19) 



(20) 



where repeated indices are summed over. Here the non-zero 
components of the matrices and are given as 

= u, mis = 2v, {uxt = ^j-v, mil = J' 



and 



IS = \j, = \j+v, mil = J' 



as a function of q in the first quadrant of the effective "large" 
Brillouin zone as well as cuts along the main symmetry di- 
rections. We compare calculations that have been performed 
correctly including the matrix elements to calculations where 
the matrix elements have been considered to be constant 
a^(/c)a^*(/c) = 1/5, as has been reported in certain ab initio 
electronic structure calculations. Here we used the five band 
fit discussed in the previous section and chose a chemical po- 
tential that corresponds to the undoped compound (x = 0). As 
can be seen from Fig. |7] the matrix elements play a very im- 
portant role by "filtering" special structures of xs{q) in mo- 
mentum space, e.g. suppressing the weight for small q. The 
matrix elements for the dx2_y2 and dsz-^-r"^ are very small 
on the Fermi surfaces (Fig. [6]) so that the actual Xs{q) is re- 
duced when these matrix elements are properly taken into ac- 
count. Neglecting the matrix elements will lead to a wrong 
result and generally to a higher value and more homogeneous 
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FIG. 8: (Color online) The RPA enhanced susceptibilities calculated 
for the electron-doped compound (x = 0.125). The interaction pa- 
rameters have been chosen as [/ = = 1.65 and J = 0. (a) and (c) 
are plots of the spin susceptibility, (b) and (d) are plots of the charge 
susceptibility. 
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FIG. 9: (Color online) The RPA enhanced spin susceptibility 
Xs{Q,(^ = 0) as a function of the interaction strength U = V for 
the electron-doped compound (x = 0.125) at two different momenta 
Q = (tt, 0) and Q* = (tt, O.IGtt). 



pound (x = 0.125). This shift can be attributed to an im- 
perfect nesting of the a and P FS sheets due to the opposite 
growth of the electron and hole sheets with doping. 

When we finally compare in Fig. [9] the change of the spin 
susceptibility as a function of U for the antiferromagnetic 
wave vector Q = (tt, 0) to the one for the incommensurate 
wave vector Q* = (tt, O.IGtt) we find a moderate and nearly 
linear increase of Xs{Q) while Xs{Q^) diverges if we ap- 
proach the critical value of U. 



distribution of the susceptibility in momentum space. Impor- 
tant features like the Q = (tt, 0)-peak, that is responsible for 
the antiferromagnetic SDW instability, can be under- or over- 
estimated. 

For the single band susceptibility the inclusion of interac- 
tions within an RPA approach is known to enhance existing 
structures in the bare susceptibility as Uxo{q) approaches 1. 
In the case of a multiorbital susceptibility with onsite intraor- 
bital repulsion U, onsite interorbital repulsion V and Hund's 
rule coupling J it is not obvious how the different structures 
in the spin and in the charge susceptibility are changed by 
the variation of these three parameters. For a simplified and 
more transparent discussion we will first study the suscepti- 
bility as a function of U and we will choose V = U and 
J = J' = 0. In Fig. [8] we compare the spin and charge sus- 
ceptibilities of the electron-doped compound with x = 0.125 
for a value ofU = 1.65 that is chosen such that the main 
peaks in the spin susceptibility are considerably enhanced. 
Here, as elsewhere in this paper, energy units are in eV. For 
this value of U we find that the charge susceptibility is more 
than one order of magnitude smaller than the spin susceptibil- 
ity. In addition the charge susceptibility has no pronounced 
structures in momentum space whereas the spin susceptibility 
shows two distinct peaks between (tt, 0) and (7r,7r). While 
in the undoped compound (see Fig. [7]) we find the main peak 
in the spin susceptibility at the antiferromagnetic wave vec- 
tor Q = (tt, 0) (corresponding to (tt, tt) in the "small" BZ) 
we observe a shift of this peak towards an incommensurate 
wave vector Q* = (7r,0.167r) in the electron-doped com- 



C. T dependence of x and (Ti T) " ^ 

The temperature dependence of the magnetic susceptibil- 
ity and NMR spin-lattice relaxation rate in the paramagnetic 
phase are immediate tests of our RPA calculation of xs (q, ^) • 
The susceptibility of the F-doped LaOFeAs 1111 material 
which we study here has been measured both by direct mag- 
netization and NMR methods^^'^^"^^ and found to increase 
quasi-linearly up to several hundred degrees K. In addition, 
Zhang et ah have pointed out that this behavior appears 
to hold for the 122 class of materials as well, and that the 
slope of the linear-T behavior appears to be roughly doping- 
independent^^. Within our theory, we can calculate the homo- 
geneous, static bare spin susceptibility by summing over spins 
and using Eq.[22l At T = it reduces to 

Xo = 2x5 = 2^7V.(0), (23) 

V 

which is of course just the Pauli susceptibility proportional 
to the total density of states at the Fermi level (A/'j,(0) is the 
single- spin density of states at the Fermi level in band v.). 
At finite temperatures one might assume that the susceptibil- 
ity of an itinerant electron system should decrease. This is 
the usual case for a single parabolic band, but in the presence 
of band structure effects the susceptibility may first increase. 
For example, in a single band system the band structure en- 
ters in a simple way as xo(T) ^ Xo(0) + [{N" {{))N{{)) - 
iV'(0)2)/7V(0)2]T2, where iV^ N" , are the derivatives of the 
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FIG. 10: (Color online) The susceptibility Xo(^) = 2xs(q = 
0,cj = 0; T) (blue) in units of eV"^Fe"^ vs. T in Kelvin for a 
doping of X = (a) and x = 0125 (b). The interaction parameters 
have been chosen to be [/ = 1.1, f = 0.25U (red) and J = J' = 0, 
t/ = y = 1.35 (a) or = y = 1.5 (b), respectively (green). 



density of states at the Fermi level. Our situation is also 
somewhat unusual as the chemical potential sits in a region 
of rapidly varying density of states, so one may expect the T 
dependence of the susceptibility to differ from the textbook 
results. 

We now plot in Fig.[TOlthe temperature dependence of both 
Xo and the RPA enhanced x for the undoped and doped cases, 
for a particular choice of interaction parameters. It is seen 
that in both the doped and undoped cases, the susceptibility 
indeed increases quasilinearly over the experimentally inter- 
esting range of 100 to 500K. This qualitative result appears 
to be relatively independent of the choice of interaction pa- 
rameters since it arises from the band structure. The RPA en- 
hancement for the cases shown is of order 30%. While this 
factor increases as one approaches the critical U at which the 
Stoner instability is reached, it cannot increase dramatically, 
because for the electronic structure of the current model, the 
instability at large q occurs first. Thus a large enhancement 
(Ref.l42h of the g = (0, 0) susceptibility does not occur in this 
model. We conclude that the observed susceptibility contains 
a significant temperature-independent interband or van Vleck 
component which is not included in the xs calculated here. 
Similar conclusions were reached by the authors of Ref.39. 

We also calculated the NMR relaxation rate l/(TiT) that is 
proportional to to the imaginary part of the local (q integrated) 
spin susceptibility 
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FIG. 11: (Color online) The NMR relaxation rate in units of 
eV~^Fe~^ as a function of temperature T for two different dopings 
and with the same set of interaction parameters and line colors as in 
Fig-Cni Panels (a) and (b) represent x = and x — 0.125, respec- 
tively, for form factor A{c\) — 1. Panels (c) and (d) represent x = 
and X — 0.125 for A(q) = cos ^cc/2 cos ^y/2. 



and (b) we take the hyperfine form factor of the ^^F nucleus 
used e.g. in Ref. 40 to be constant, since the F is located di- 
rectly above the Fe in the LOFFA material. In panels (c) and 
(d), we have modelled the form factor of the As nucleus as 
74(q) = cos ^2^/2 cos g'y/2, which suppresses the (tt, 0) fluctu- 
ations, since the As atom is in the center of a square containing 
four Fe. 

Figures [TOl and [TT] are intended to show some of the de- 
pendence of the susceptibility and T^^ on interaction param- 
eters for a few cases. The trends are qualitatively similar to 



experiment: 



,38,39.40 



, although the increase with increasing tem- 



perature for both quantities found here is somewhat weaker 
than observed. We have not attempted to fit experiments in 
detail. 



where A(ci) is a geometrical structure factor and cjl is the 
NMR frequency. In Fig[TT]we show l/(TiT) as a function 
of temperature for the same set of parameters that we have 
used for the q — ^ susceptibility. Here we find for the in- 
teracting system an initial decrease of the relaxation rate as 
a function of temperature that is most significant for the un- 
doped compound with finite values of J. In the latter case we 
find a very distinct upturn around 300 K followed by a slight 
and nearly linear increase. For the other cases we find a sim- 
ilar tendency, although not as pronounced. In panels [TT] (a) 



V. THE EFFECTIVE INTERACTION AND THE GAP 
FUNCTION 



Assuming that the pairing interaction responsible for the 
occurrence of superconductivity in the iron pnictides arises 
from the exchange of spin and charge fluctuations, we can 
calculate the pairing vertex using the fluctuation exchange ap- 
proximation^^. For the multi-orbital case^^, the singlet pairing 
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vertex is given by 



tq 



(25) 



The Xi^^^ term describes the spin fluctuation contribution 
and the Xo^^"^ t^™ the orbital (charge) fluctuation contribu- 
tion. Fig. [12] shows the Matsubara frequency dependence of 
these two terms for a momentum transfer (tt, 0) and a typical 
interaction strength. Here one sees the that the spin fluctua- 
tion contribution is dominant and falls off on a frequency scale 
which is small compared with the bandwidth. Thus while the 
gap equation depends upon the kernel Imr^^(/c, k' ^uo)), the 
important k and k' values are restricted by this frequency cut- 
off to remain near the Fermi surfaces. Then, just as for the 
electron-phonon case, the strength of the pairing interaction is 
characterized by a frequency integral of this kernel weighted 
by uj~^. Making use of the Kramers-Kronig relation 



7TUJ 



= Re[rl^{k,k\u; = 0)] (26) 



we can proceed further by considering only the real part of the 
CO = pairing interaction. If we now confine our considera- 
tions to the vicinity of the Fermi surfaces we can determine 
the scattering of a Cooper pair from the state {k^—k) on the 
Fermi surface Ci to the state {k' ^—k') on the Fermi surface 
Cj from the projected interaction vertex 

r,,(fc,fc') = ^a*f(-fc)<;*(fc)Rere(fc,fc',0)] 



stpq 



xalXk')ali-k') 



(27) 



where the momenta k and k' are restricted to the different 
Fermi surface sheets with k e Ci and k^ G Cj. If we de- 
compose the superconducting gap into an amplitude A and a 
normalized symmetry function g{k) we can define a dimen- 
sionless pairing strength functional^"^ as 



A[ff(fc)] 



Here Tij is only the symmetric part 



2 i^'i^ji^^ k ) -\- r^j(fe, k ) 



(28) 



(29) 



of the full interaction, which gives identical results within 
the spin singlet subspace. The Fermi velocity is defined to 
be VF{k) = \\/kEiy{k)\ for k on the given Fermi surface. 
From the stationary condition we find the following eigen- 
value problem 



E 



dk' 



27r 2i:vF(k') 



-r,,(fc,fc')5a(fc') = (30) 







spin + orbital — • — 
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FIG. 12: (Color online) The energy dependence of the spin and or- 
bital part of the effective pairing interaction r(/c, /c^ cjn, cj^)- 



The kernel Vij(k^k') is evaluated at temperatures well be- 
low the characteristic temperature at which the spin fluctu- 
ation spectrum has formed. In this temperature region, the 
interaction is independent of temperature. From the above 
eigenvalue problem we will determine the leading eigenfunc- 
tion ga{k) and eigenvalue Aq, for a given interaction vertex 
Vij {k^k'). The largest eigenvalue will lead to the highest tran- 
sition temperature and its eigenfunction determines the sym- 
metry of the gap. The next leading eigenvalues and eigen- 
functions further characterize the pairing interaction and can 
indicate the structure of possible collective states. With the 
knowledge of the pairing function ga{k) we can also deter- 
mine the individual contributions of the different intra- and 
interorbital scattering processes to the total pairing strength 



VI. RESULTS FOR THE PAIRING STRENGTH A AND 
THE GAP FUNCTION g{k) 

A. Eigenvalue problem 

In this section, we present and discuss results for the pairing 
strength A and the symmetry function g{k) obtained within 
our weak-coupling approach. 

Solving the eigenvalue problem given in Eq. [30l we find 
a set of eigenf unctions ga{k) defining the k dependence of 
the gap on the FS sheets together with a set of corresponding 
eigenvalues Aq, denoting the dimensionless pairing strengths 
associated with a given ga{k). We first have to classify the 
different eigenfunctions according to basic symmetry opera- 
tions. In the following we will speak of an s wave state if 

9{~kxjky) = g{kxj —ky) = g[kx^ ky) 
giky.kx) = g{kx,ky) (31) 

while we have a (ia.2_^2 wave state if 

9{~kxjky) = g{kxj —ky) = g[kxj ky) 
g{ky,kx) = -g{kx,ky) (32) 

Since there are two different dx'^-y^ wave states among the 
first few eigenfunctions, we have to distinguish them further. 
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e.g. by comparing the sign of gu{k) for u = ai and u = a2 
in the same direction in momentum space. Here we will label 
the state that changes sign between ai and 0^2 with dj.2_y2 (1), 
the one without sign change with dx2_y2{2). Furthermore we 
can distinguish didxy wave state that is given by 



gikyi kx) = gikx") ky) 



and a g wave state with 

g{~kx-,ky) = 
g{ky^ kx) = 



g{kx-, ky) 
~g{kxj ky) 



g{kx^ ky) 



g{kx-, ky) 



(33) 



(34) 



B. Results: J = J' = 

First we consider a case where the Hund's rule coupling and 
the pair hopping energy are negligible compared to the intra- 
and interorbital Coulomb interactions, so we set J = J' = 
and V = U.ln Fig. [13] (a) we show the pairing strength eigen- 
values \a for the four leading eigenvalues as a function of U 
for the electron-doped compound (x = 0.125). Approaching 
the critical value of U where the eigenvalues start to diverge 
we find a clear separation of the two leading eigenvalues from 
the next two eigenvalues. However the two leading eigenval- 
ues, corresponding to the s and the dx2-y2 symmetry remain 
very similar in size and we find a crossover from the dx2-y2 
to the s symmetry around U = 1.65. We also show the s wave 
(d) and the dx2-y2 wave (e) pairing functions on the four FS 
sheets foxU = 1.73 close to its critical valued. The extended 
s wave state (which we have labelled s) is characterized by i) 
a sign change on the 13 FS sheets with nodal points displaced 
from the generic (0, 7r)-(7r, 0) direction that would result from 
a pure cos kx^cos ky state, and ii) a change of average sign on 
the [5 sheet relative to that on the a sheets. On the a FS sheets 
we find a nodeless but anisotropic gap distribution with higher 
weight on the small ai compared to the larger 0^2 sheet. The 
dx2-y2 State features an anisotropic gap distribution on the f3 
FS sheets and a rather conventional d wave gap distribution 
on the a FS surfaces, with a sign change between the ai and 
a2 sheet. In Fig. [14] we show the corresponding results for 
the undoped (x = 0) compound. Here we find similar results, 
although the eigenvalues close to the instability are better sep- 
arated and the crossover from dx2-y2 to s appears already for 
a rather small value of /7. If we compare the s wave symmetry 
function for x = close to the instability (see Fig. [14] (d)) to 
the corresponding symmetry function for the electron-doped 
compound, we find that the nodal points on the (3 FS have 
moved even closer to the tips of the sheets and the negative 
dip between them has decreased considerably. For the dx2-y2 
symmetry (see Fig. [14] (e) we find that the weight on the 0^2 
FS sheet has nearly vanished. 

In Figs. [13] and [14] as well as in the following Figs. [T6] and 
[T5] we compare in panels (b) and (c) the contributions of the 
different intra- and interband processes. Here \oicx sums up 
all contributions of \ij resulting from scattering within each a 
sheets and in between the two a sheets as X^a = Xli 7=12 
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FIG. 13: (Color online) The eigenvalues and eigenfunctions for the 
electron-doped compound {x — 0.125) ioxU — V and J — J' — 
The four largest eigenvalues as a function of IJ (a) and the different 
inter- and intraband contributions to the eigenvalues A for the two 
symmetries with largest eigenvalues, extended s (b) and d^2_y2 (c) 
wave. Color coded plot of the extended s wave (d) and the d^2_y2 
wave (e) pairing functions along the different Fermi surface sheets, 
calculated close to the instability {U — V — 1.73). 



where i = 1 refers to the ai, i = 2, to the 0^2, i = 3 to the 
(3i and z = 4 to the (32 Fermi surface. A^^^^^ = A33 + A44 
denotes the intraband contributions of the $ FS sheets, while 
X'tnter _ _^ ^j^^ Corresponding interband contribu- 

tion. Finally Xap sums up all the remaining contributions, re- 
sulting from scattering between the a and (3 Fermi surfaces. It 
is obvious that in the case of J = J' = for both of the main 
pairing symmetries the intersheet Xap contribution is respon- 
sible for the rapid increase of the pairing strength A with U, 
reflecting the nesting between the a and (3 FS sheets with nest- 
ing vector Q*. All other contributions are small and mainly 
negative. For finite values of J and J\ which we will next 
discuss we find that for the dx2-y2 pairing symmetry the in- 
terorbital contributions between the two (3 sheets become also 
important and for the electron-doped compound even domi- 
nant. 
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FIG. 14: (Color online) The same as in Fig.[T3]but for the undoped 
compound (x — 0) with V — U and J — J' — ^. Here the two 
pairing functions shown in (d,e) correspond to the two leading eigen- 
values and are calculated for [/ = 1.54. 



FIG. 15: (Color online) The same as in Fig.[T3]for the undoped com- 
pound (x = 0) with J' = J/2 = [7/4 and = [/ - 5/4 J. Again 
the extended s wave (d) and the dx2_y2 wave (e) pairing functions 
are calculated close to the instability with U = 1.25. 



C. Results: J, > 

Now we consider a case of finite J and J^ Here we will 
choose = J/2 and we will fix F = U - 3/4J - J'. 
These choices are consistent with the generation of intrasite 
couplings from a Hubbard argumentf^, as mentioned above, 
as well as with the range of values for J/U found by Anisi- 
mov et ah in an ab initio calculation^"^. In Fig.[T5](a) we show 
for the undoped compound {x = 0) and for = J/2 = [//4 
the same four eigenvalues as a function of U as in the pre- 
vious figures. Here we find again that the extended s wave 
state (d) is the most stable pairing configuration followed by 
\hQdx2_y2 wave state (e). For iht dx2_y2 state the main con- 
tributions to the pairing function are along the (3 FS sheets, 
while the a FS sheets contribute less significantly. Since this 
state has no nodes along the 13 sheets, it can be considered 
as a mainly nodeless state, in contrast to the states found for 
J = J' = 0. The same is true for the electron-doped case 
with X = 0.125 (Fig. [16]). Here the dx2_y2 state is the most 



stable state followed by the extended s wave state. Close to 
the instability the nodeless P sheets contribute dominantly to 
the pairing. 

Although the pairing is constrained to the close vicinity of 
the Fermi surfaces, one can try to find an approximation of 
the gu{k) on the different FS sheets that extends further into 
the Brillouin zone. Here we want to restrict ourselves to the 
smallest number of harmonics necessary to find a reasonable 
fit of the pairing function on the different FS sheets. For the s 
wave symmetry we can write 

gu{k) = 2 [cos kx -\- cos ky -\- 2wiy^xy ^oskx cos ky 

-\-2wiy^4x4:y COS Akx COS Aky] (35) 

where we find for [7 = 1.5 the following parameters: Aq,^ = 
0.051, Wai,xy = —0.35, and Wai,4x4y = —1.35 for the ai 
sheet, = 0.02, Wa2,xy = —0.1, and Wa2,4x4y = 1 for the 
a2 sheet, and Ap = 0.15, wp^xy = 0.17, and Wi3^4x4y = —0.1 
for the P sheets. In the same way, we can find an approxima- 
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FIG. 16: (Color online) The same as in Fig.[T3]for the electron-doped 
compound (x 0.125) with f J/2 U/4 and F = [/ - 5/4J. 
Here the dx2_y2 wave (d) is more stable than the extended s wave 
(e) pairing function and the corresponding eigenvalues are clearly 
separated. Both are calculated close to the instability with U = 1.28. 



tion to the d wave symmetry function as 

gu{k) = 2Ajy [cos kx — cos ky + Wjy 2x(cos 2kx — cos 2ky)] 

(36) 

with = — 1.2 and Wai,2x = for the ai sheet, A^^ 
0.1 2 and i(;q;2,2x = 0, for the 0^2 sheet, and Ap = —0.018 and 
'^/3,2cc = —1.12 for the sheets. The results of the fitting are 
shown in Figs. [TTl and [TSl 



D. Role of nesting 

To study the influence of nesting on the spin susceptibil- 
ity, especially on the (tt, 0) peak, we slightly modify the hop- 
ping parameters, creating a toy model with perfectly nested 
FS sheets. This means that we try to find an approximation 
where the 0^2 and the f3 FS sheets are of approximately the 
same size and shape. The band structure used within this toy 
model is very similar to the band structure found formerly by 




0.2 0.4 0.6 0.8 1 

FS 



0.2 0.4 0.6 0.8 1 
FS 




0.2 0.4 0.6 0.8 1 
FS 



0.2 0.4 0.6 0.8 1 
FS 



FIG. 17: (Color online) The s wave symmetry function g{k) on the 
different FS sheets forU = 1.5 and x = 0. Here we compare a 
fit of gu{k) (blue line) to the actually calculated values of g{k) (red 
points). 
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FIG. 18: (Color online) The d wave symmetry function g(k) on the 
different FS sheets forU = 1.5 and x = 0. Here we compare a 
fit of gu{k) (blue line) to the actually calculated values of g{k) (red 
points). 



12 




FIG. 19: (Color online) The FS sheets found within a circular ap- 
proximation showing nearly perfect nesting between the 0:2 and the 
P FS sheets with nesting vectors (tt, 0), or (0, tt), respectively. Here 
we use the same color convention as in Fig.|5] 
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FIG. 20: (Color online) The susceptibility xo {q) for a toy model with 
spherical FS sheets. 



a fit to the DFT bands, to assure that all basic properties of the 
model are reasonable, i.e. that the matrix elements are simi- 
lar to the correct matrix elements found from the DFT fit. In 
Fig. [19] we show the near circular FS sheets of this model. 

The bare spin susceptibility Xo(^) for this model in the un- 
doped case is exhibited in Fig. [20l and may be compared to 
the results of the original model in Fig. [71 We see that the 
(tt, 0) peak remains but is not particularly enhanced by the en- 
hanced nesting, implying that the original model was already 
quite close to a nesting condition. On the other hand, some 
subtle differences which are quite interesting appear when we 
examine the pairing functions. If we consider an even more 
Gedanken-typc model where the two a sheets have degener- 
ate radius a, which is also the same as the two (3 sheets, we 
see that the Fermi surface in the 1st effective Brillouin zone is 
invariant under a translation by (tt, 0). Under this transforma- 
tion, simple extended-s and (1^2 _y2 functions cos kx + cos ky 
and cos kx — cos ky map into one another identically. We 
therefore expect that the pairing eigenvalues for s and dx'2-y2 
will become degenerate. The toy model band structure con- 
sidered here and shown in Fig. [19] is nearly the same as the 
Gedanken model, but has two slightly nondegenerate a sheets. 
We see nonetheless in Fig.[2T](b,c) and (d,e) that the two com- 
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FIG. 21: (Color online) The leading eigenvalues A for the toy model 
with spherical FS sheets (a). The two leading pairing functions g(}z) 
for the toy model with spherical FS sheets for [/ = V — 1.6 (b,c) 
and for = y = 1.66 (d,e). 



peting states are indeed nearly exactly degenerate. Fig. [21] (a) 
exhibits the dependence of the leading eigenvalues on IJ . 



VII. THE SPATIAL AND ORBITAL STRUCTURE 

The gap function gv(k) contains information on the spatial 
and orbital structure of the pairs. In the previous section, we 
have calculated gv(k) and found that its behavior on the Fermi 
surfaces can be fit by low order harmonics for both the s and 
d-wave eigenstates. Here we use these results to determine 
simple pictures of the internal pair structure. A gap operator 
can be written as 

^=^E5>^(^)^-TW7.i(-fc) (37) 



v.k 



with jua{k) the destruction operator for an electron in the z/^^ 
band with wave vector k and spin a. Using the band-orbital 
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matrix elements to relate the band operator to the orbital op- 
erator 

lua{k) = Y,a]^{k)dna{k) (38) 

n 

we have 

A = ^ ^ anmih - k)dn]{h)dmi{k) (39) 

with 

= ^ ^5.(fc)a;!(fc)a-(-fc)e-*' (40) 

ku 

where I = I1 — I2. The amplitude anm (0 describes the internal 
spatial and orbital structure of the pair. 

Using the harmonic approximations for gu{k) for the s and 
dj;2_y2 gaps discussed in the previous section, we have cal- 
culated anm{l)- In principle the k sum should be cut-off 
when k — kp exceeds the inverse coherence length. Here we 
have used a Gaussian cutoff of the k sum which decays when 
|k — ki?| exceeds 27r/^o, where ^0 is taken to be of order 3 
times the Fe-Fe spacing a. This provides a local picture of the 
internal orbital structure of a pair. This basic structure contin- 
ues out to a radius set by the coherence length ^0 • 

We find that the off-diagonal amplitudes anm are negligible 
compared to the diagonal ones, and that the orbitals that con- 
tribute are the dxz, dyz, and dxy orbitals which have weight 
near the Fermi surfaces. The amplitudes ann{l) for the s and 
the dx2-y2 gaps are shown in Figs. [22land[23l respectively. 
For the s case, one sees that the internal structure of a pair 
consists of a superposition of {xz ]^xz [) singlets which are 
formed between a central site and sites displaced by an odd 
number of lattice spacings in the ^/-direction, {yz ]^yz [) sin- 
glets between the central site and odd numbered sites in the x- 
direction, and weaker {xy ]^xy [) singlets with a more intri- 
cate s-wave arrangement. The internal structure of the d^'^-y^ 
state consists of a similar superposition but with a negative 
phase difference between the {xz ]^xz [) and {yz ]^yz [) 
singlets. In addition, there is a significant dx2-y2 contribu- 
tion from the {xy ]^xy [) singlets confined primarily to the 
nearest neighbor sites. 

VIII. SUMMARY AND CONCLUSIONS 

We have studied a weak-coupling spin and charge fluctu- 
ation model of the pairing interaction using a tight-binding 
parametrization of the electronic band structure obtained from 
DFT calculations for the Fe-pnictides. We initially reviewed 
criticism of the 2-orbital models used in early studies of this 
kind, and compared results with more accurate 5 -orbital mod- 
els. Within the 5 -orbital framework, we have calculated the 
multiorbital susceptibility. We calculated the static homoge- 
neous spin susceptibility within this framework, and showed 
that it qualitatively agrees with experimental measurements of 
the T dependence, but that an interband susceptibility com- 
ponent was necessary to understand the magnitude of the 
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FIG. 22: (Color online) The spatial and orbital pair structure for the 
s gap calculated from the harmonic approximation of g{k) corre- 
sponding ioU — 1.5 with a cut-off length of Ac = 3.3dFe-Fe. 



measured low-temperature limiting susceptibility. The T- 
dependence of both x(^) {TiT)~^ is qualitatively similar 
to experiments, and may be consistent with them for some val- 
ues of the parameters. A detailed fit was not attempted here. 

We then constructed the pairing vertex from the general- 
ized RPA susceptibilities, and calculated the pairing eigenval- 
ues for dopings corresponding to undoped and electron-doped 
materials, for a variety of interaction parameters correspond- 
ing to intra- and inter-orbital Coulomb matrix elements and 
Hund's rule couplings. We found that in the parameter region 
where there is a significant coupling strength, the pairing in- 
teraction for our model of the Fe-pnictides arises from the ex- 
change of spin fluctuations. These give rise to both intra- and 
inter-Fermi surface scattering processes. In a number of cases 
the dominant pairing contribution came from particle-particle 
scattering processes from the hole Fermi surface around the 
F point to the electron Fermi surfaces around the (tt, 0) and 
(0, tt) points in the unfolded Brillouin zone. These scatterings 
involve momentum transfer of order (tt, 0). We find that there 
are two leading pairing channels, one with "sign-changing" s- 
wave {Aig) symmetry and one with dx2-y2 {Big) symmetry, 
and that the gap functions corresponding to both have nodes 
on the Fermi surface. For values of the interaction parameters 
corresponding to significant pairing strength, these two eigen- 
values can become very close. The s-wave gap that exhibits 
a sign change between its average values on the electron and 
hole Fermi surfaces nevertheless also displays nodes on the 
electron sheet, whereas the d-wave state has its nodes on the 
hole sheet. The s-wave gap nodes are not required by symme- 
try, and could be absent for a different choice of parameters, 
but appear to be robust within the manifold of apparently re- 
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FIG. 23: (Color online) The spatial and orbital pair structure for the 
dx2_y2 gap calculated from the harmonic approximation of g{k) cor- 
responding to [/ = 1.5 with a cut-off length of Ac = 3.3dFe-Fe- 



alistic interaction parameters we have studied. 

An obvious question which arises in such studies is the sen- 
sitivity of results to the particular choice of band structure. 
We have relied upon the Fermi surface for the paramagnetic 
system determined in the Generalized Gradient Approxima- 
tion (GGA) calculations of Cao et alA, but other band calcu- 
lations find subtle differences in Fermi surfaces for the La- 
1111 material. We have attempted to address this question by 
considering variations of our effective 5 -band model hopping 
parameters chosen originally to fit the Cao et ah band struc- 
ture. In particular, we considered variations which led to a 
perfect nesting of the outer a and 13 Fermi surface sheets, to 
try to maximize the (tt, 0) contribution to the pairing vertex. 
We found that our results for the structure of the susceptibility 
and pairing vertices changed in fact very little. On the other 
hand, the subtle changes led to a nearly exact degeneracy of 
the extended-s and dj;2_y2 eigenvalues of the linearized pair- 
ing problem. We argued that this degeneracy becomes exact 
for a situation where all sheets have the same radius, and pro- 
pose this as a simple explanation of our finding that these two 
pairing channels appear to compete very strongly with one an- 
other, even in the realistic cases. 

Our results appear to be qualitatively different from several 
^Qj.]^^i5ji8j2J;6j47 f2j^(^^ on heuristic grounds or within 

microscopic calculations, that the Aig state is closer to the 
form cos kx x cos ky (in the unfolded zone), in contrast to our 
result, which is closer to the form cos kx + cos ky in the sense 
described above. Although both states have identical Aig 
symmetry, the difference between them is important given the 
structure of the Fermi surface; for example, the cos kx x cos ky 
form does not have nodes on either of the p sheets shown in 



Fig.m We believe that the strong pair scattering between the 
f3i and (32 sheets found in our calculations, and neglected in 
some simpler models, is crucial to stabilize the gap structure 
similar to cos kx -\- cos ky. 

We further studied the spatial structure of the states corre- 
sponding to the leading pairing eigenvalues. In each case, the 
internal structure consists of a superposition of singlets made 
up from electrons which occupy the same dxz, dyz, or dxy or- 
bitals on different sites. The distribution and relative signs of 
these superpositions are illustrated in Figs.[22]and[23l The fact 
that the s-wave and d-wave solutions have very similar cou- 
pling constants opens the possibility that different members 
of the Fe-pnictide superconducting family may have different 
gap symmetries. Furthermore, it suggests that these systems 
may have low-lying collective particle-particle modes which 
would be 5-wave like in a d-wave superconductor, and d-wave 
like in an s-wave superconductor"^^. 

There have been several earlier weak-coupling calculations 
of pairing in the Fe-pnictide superconductors. An early sug- 
gestion of an extended 5-wave sign reversed gap made by 
Mazin et al^^ was based upon the Fermi surface structure 
found in a DFT calculation. There it was proposed that (tt, 0) 
antiferromagnetic spin fluctuations could lead to a supercon- 
ducting state with an isotropic s-wave gap which had opposite 
signs on the electron and hole pockets. 

Using a 5 -orbital parameterization of the DFT bandstruc- 
ture, Kuroki et al}^ have carried out a calculation in which an 
RPA form for the spin and orbital contributions to the inter- 
action was used along with bare single particle Green's func- 
tions to construct a linearized gap equation. In a parameter 
range similar to ours, they found that the leading pairing in- 
stability had 5-wave symmetry with nodes on the electron-like 
Fermi surface. They also noted that the next leading channel 
had dx'2-y2 symmetry and discussed conditions under which 
this could become the leading pairing channel. The behavior 
of their s-wave gap differs from what we find by a phase fac- 
tor of -1 on the electron Fermi surfaces. That is, if the gap 
on the inner hole Fermi surface around the T point is taken 
as positive, then Kuroki et al find that the gap on the elec- 
tron Fermi surface at (tt, 0) reaches its largest negative value 
at the point closest to the origin F. With the same convention, 
we find that the gap on the electron Fermi surface takes on its 
largest positive value at this point. This difference may reflect 
a difference in the nesting of the hole and electron Fermi sur- 
faces due to fits to slightly different band structures or from 
the choice of interaction parameters'^^. 

Wang et al}^ have studied the pairing problem using a 
5 -orbital effective band structure together with the func- 
tional renormalization group approach. Using the same tight- 
binding parameterization as Kuroki et al , they also find that 
the leading pairing channel has 5-wave symmetry and that the 
next leading channel has dx'2-y2 symmetry. For their interac- 
tion parameters, they find that there are no nodes on the Fermi 
surface, but there is a significant variation in the magnitude of 
the gap. With the sign convention where the gap on the inner 
hole Fermi surface is positive, their gap function reaches its 
smallest negative value at the point on the electron Fermi sur- 
face which is closest to V. They note that for other parameter 
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choices, this 5- wave gap function on the electron Fermi sur- 
face may have nodes. This would lead to a gap which becomes 
positive on a region of the electron Fermi surface closest to the 
r point, in agreement with what we have found. 

Both of these studies also conclude that the pairing interac- 
tion arises from spin fluctuation scattering with the dominant 
contribution associated with Q ~ (tt, 0) scattering of a pair 
from the inner Fermi hole surface around the F point to the 
electron Fermi surface around (tt, 0) and (0,7r). They also 
find that the magnitude of the gap on the central hole Fermi 
surface is smaller than the maximum magnitudes of the gap 
on the inner hole and the electron Fermi surfaces. Thus the 
five-orbital weak coupling calculations all find singlet pairing 
in the extended s-wave channel with a nearby dx2_y2 wave 
state. The question of whether there are gap nodes for the s- 
wave, and indeed whether or not the s-wave state is ultimately 
stable with respect to the d-wave state, would both appear to 
depend on parameters. This sensitivity appears to us to be a 
natural consequence of the importance of several orbitals near 
the Fermi surface. In contrast to the cuprates, where calcu- 
lations of this kind show a clear dj;2_y2 state well separated 
from other pairing states, it seems possible here that variations 
in band structure or interaction parameters found in different 
materials might possibly lead to different symmetry ground 
states. 

We close by pointing out that at this writing, experiments 
have not conclusively answered either the question of order 
parameter symmetry or even whether gap nodes are present. 
Early indications from specific heat^, Andreev point con- 
tact spectroscopy^^, and NMR^ on the 1111 materials sug- 
gested nodes in the superconducting order parameter. On 
the other hand, some penetration depth experiments^^ and 
ARPES experiments^^'^^ appear to find nearly isotropic gaps 
in the 122 materials. Our calculations find nearly degenerate 
s and dx2_y2 states at Tc, both of which have nodes on the 
Fermi surface. If the experiments indicating a lack of nodes 
are correct, there are two ways to imagine reconciling our con- 
clusions with experiments on the 122 materials. The first is to 
consider the effects of disorder present at fairly high levels in 
current crystals, which should average the order parameters to 
a finite quasiisotropic value in the s case. The second possi- 
bility follows from the observation that if the dj.2_y2 state is 
the leading eigenvalue at Tc, the thermodynamic ground state 
may be unstable towards an admixture of the s state^^ . A simi- 
lar phenomenon occurs in the phase diagram of ordinary s and 



d symmetries in one-band superconductors with these pairing 
channels competing^^. Separating these possibilities may be 
difficult, and will probably require phase sensitive probes as 
were discussed in the cuprate context^. 
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APPENDIX: MODEL PARAMETERS FOR THE 5 BAND FIT 

In the following, the model parameters for the 5 band tight 
binding fit of the DFT band structure by Cao et ah are listed. 
Including hopping up to fifth nearest neighbors on an effective 
Fe-Fe lattice we find intraorbital kinetic energy terms ^rnm {k) 
and interorbital kinetic terms ^mn(^) that respect the basic 
symmetry requirements imposed by the point group of the 
crystal. The intra- and interband hopping parameters are given 
in Tab. Il and Tab. [Ill respectively. The onsite energies (in 
eV) ei = 0.13, 63 = -0.22, 64 = 0.3, and 65 = -0.211 
where i = 1 corresponds to the dxz, i = 2io the dyz, i = 3 
to the dx2_y2, z = 4 to the dxy, and i = 5 to the d^z'^_r'2 or- 
bital, and the band structure is listed below. For the toy model 
with nearly spherical Fermi surface sheets we have changed 
tl\ from 0.28 to 0.3 and tf from -0.354 to -0.409. This is 
sufficient to produce nearly degenerate and spherical 0^2 and 
[5 FS sheets. For the Coulomb interaction strengths we have 
used, the bare ratios U /t]^y and U /t]^ are of order 5, compa- 
rable to values ofU /t 'm the cuprates. However, we believe 
that for the Fe-pnictide system, the interaction is "spread out" 
across multiple bands leading to a subtantially less strongly 
correlated system^. 
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^11/22 = 24/^ COS kx + 2t^/^ cos ky + 4^^ cos /c^, cos ky ± 24i(cOS - cos 2ky) + ^t]^xy/xyy COS COS /c, 

COS 2ky COS fe^ + 4^^^^^ cos(2fe^) cos(2fe^) 
^33 = 24^(cos kx + COS ky) + At^^y cos A;^, cos A;^^ + 24^(cos 2kx + cos 2ky) 
^44 = 2t^^(cos /Ca^ + COS ky) + 4t^^ COS kx COS /cjy + 2t^^(cos 2kx + COS 2/c^) 

+4t^^^(cOS 2kx cos /C^y + cos 2ky cos /C^^) + ^t^tyy COS 2/Ca^ COS 2/^2; 

^55 = 2tl^{coskx + cosA;^^) + 24^(cos2A:^ + cos2ky) 

+4t^^^(cOS 2kx cos ky + cos 2ky cos /C^^) + ^t^xxyy COS 2/Ca^ COS 2/^2; 

^12 = -44^ sin /ca, sin ky - At]^^y{sm 2kx sin ky + sin 2/c^ sin kx) - At]^xyy si^i 2/ca; sin 2ky 

^"i3/23 = ±2z4^ sin ky/x ± 4z4^ sin ky/x cos /c^/^ T 4z4^^(sin 2/c^/^ cos /c^/^ - cos 2kx/y sin /c^/^,) 

^14/24 = 2^4^ sin /Ca^/^ + 4^4^ COS kyjx sin /c^,/^ + ^i^t^y sin 2/Ca;/^ cos /c^/g, 

^15/25 = 2^4^ sin /c^/a, - 4^4^ sin ky/x cos kx/y - ^itHyy sin 2/c^/^ cos 2kx/y 

^34 = 44^2/ (^^^ "^^y ~ ^2/) 

^35 = 24^ (cos kx — cos /c^) + 4:tll.y{cos 2kx cos /c^y — cos 2ky cos /c^^) 

^45 = 44^ sin kx sin ky + ^t^xyy sin 2/ca^ sin 2/c^ 



j_mm 


i — X 


i = y 


i = xy 


i = XX 




i = xyy 


i — xxyy 


m = 1 


-0.14 


-0.4 


0.28 


0.02 


-0.035 


0.005 


0.035 


m = 3 


0.35 




-0.105 


-0.02 








m = 4 


0.23 




0.15 


-0.03 


-0.03 




-0.03 


m — 5 


-0.1 






-0.04 


0.02 




-0.01 



TABLE I: The intraorbital hopping parameters used for the DFT fit 
of the 5 orbital model. 



j.mn 


i — X 


i — xy 


i — xxy 


i — xxyy 


mn = 12 




0.05 


-0.015 


0.035 


mn = 13 


-0.354 


0.099 


0.021 




mn — 14 


0.339 


0.014 


0.028 




mn — 15 


-0.198 


-0.085 




-0.014 


mn = 34 






-0.01 




mn = 35 


-0.3 




-0.02 




mn — 45 




-0.15 




0.01 



TABLE IL The interorbital hopping parameters used for the DFT fit 
of the 5 orbital model. 
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